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Abstract 

Electron-electron interactions can induce Fermi surface deformations which break 
the point-group symmetry of the lattice structure of the system. In the vicinity of 
such a " Pomeranchuk instability" the Fermi surface is easily deformed by anisotropic 
perturbations, and exhibits enhanced collective fluctuations. We show that critical 
Fermi surface fluctuations near a d-wave Pomeranchuk instability in two dimensions 
lead to large anisotropic decay rates for single-particle excitations, which destroy 
Fermi liquid behavior over the whole surface except at the Brillouin zone diagonal. 
PACS: 71.10.Fd, 74.20.Mn 

Recently a number of authors pointed out that the Coulomb interaction between 
electrons in a metal can lead to Fermi surface deformations which break the orientational 
symmetry of the system. Referring to a stability criterion for normal Fermi liquids by 
Pomeranchuk, we use the term "Pomeranchuk instability" (PI) for such symmetry- 
breaking Fermi surface shifts. Deformations of d^^-type, for which the Fermi surface 
expands along the k x -axis and shrinks along the fc y -axis (or vice versa) , have been found for 
various model Hamiltonians on a two-dimensional square lattice: t-J,|2. Hubbard, |31S E] 
and extended Hubbard model. jHj For a Fermi level close to the van Hove singularity the PI 



occurs already in the weak-coupling regime, where perturbative methods may be applied. 
Symmetry-broken Fermi surfaces in fully isotropic (not lattice) two- and three-dimensional 
Fermi liquids have also been considered. [7] From a pure symmetry-group point of view a 
PI leads to a "nematic" electron liquid as defined by Kivelson et al. [Hj in their discussion 
of a possible correspondence between electron states in doped Mott insulators and liquid 
crystal phases. 

Symmetry-breaking Fermi surface deformations generally compete with other insta- 
bilities, but may also coexist with other types of symmetry-breaking order. For example, 
a superconducting state with a d-wave deformed Fermi surface [9 J is stabilized in the two- 
dimensional Hubbard model with a sizable next-to-nearest neighbor hopping amplitude 
and an electron density near van Hove filling, at least at weak coupling. [H] Supercon- 
ducting nematic states have also been considered as one among several possibilities in a 
general symmetry classification by Vojta et al. JU] I n the following we will however focus 
on symmetry-breaking Fermi surface deformations in an otherwise normal state. 

For electrons on a lattice, the PI breaks only a discrete symmetry, the point-group 
symmetry of the lattice. Hence, no Goldstone mode exists and symmetry-broken states 
can exist also at finite temperature in d > 2 dimensions. Order parameter fluctuations 
thus suppress the PI much less than competing instabilities toward states which would 
break a continuous symmetry. In the ground state a PI can be driven as a function of 
the electron density or other control parameters, leading to a quantum critical point at 
the transition between the symmetric and symmetry-broken states, provided that no first 
order transition occurs. 

An interacting electron system in the vicinity of a PI is characterized by a "soft" 
Fermi surface, which can be deformed very easily, that is at low energy cost. In this 
article we show that strong dynamical fluctuations of such a soft Fermi surface affect 
physical properties of the system very strongly, and lead in particular to a very fast decay 
of single-particle excitations. The anisotropy of the Fermi surface fluctuations in a lattice 
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system leads to a pronounced anisotropy of the single-particle decay rate. For d x 2_ y 2-wme 
fluctuations the decay is maximal near the k x - and A^-axes and minimal near the diagonal 
of the Brillouin zone. 

The "softness" of a Fermi surface can be quantified by the Fermi surface susceptibility 

m 

which measures the Fermi surface shifts 5s kF for small momentum dependent shifts of 
the chemical potential <5/i k ' F at points on the Fermi surface. Close to a PI one of the 
eigenvalues of the matrix K ki , k ' F diverges, and the corresponding eigenvector describes the 
shape of the incipient Fermi surface deformation. 

It is intuitively plausible that critical Fermi surface fluctuations will strongly affect 
single-particle excitations and as a consequence the low energy properties of the system. 
To explore the resulting physics, we define and analyze a phenomenological lattice model 
with an effective interaction chosen such that a PI occurs, but no other instabilities. The 
model Hamiltonian reads 

H = Ck n ^ + ^7 Yl /kk'(q) n k(q)ri k /(-q) , (2) 

k,<r k,k',q 

where e k is a single-particle dispersion, n k (q) = J2 a c !c-q/2 o- c k+q/2,o-) and V the volume 
of the system. Since the PI is driven by interactions with small momentum transfers, 
that is forward scattering, we choose a coupling function / kk '(q) which contributes only 
for relatively small momenta q. This suppresses other instabilities such as superconduc- 
tivity or density waves. We emphasize that this model is adequate only if Fermi surface 
fluctuations are the dominant fluctuations in the system. Otherwise it would have to 
be supplemented by other interactions with large momentum transfers. The interplay of 
Fermi surface and other fluctuations opens a wide field for investigations in the future. 
For an analytical treatment we assume that the momentum dependence of the coupling 
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function in (2) is separable, that is 



/kk'(q) =s(q) d k d k , . 



(3) 



Although the above model can be defined in any dimension, we now focus on the par- 
ticularly interesting case of a two-dimensional system on a square lattice. To gener- 
ate a PI with d x 2_ y 2 symmetry, the form factors d^ must have that symmetry, such as 
dk = cosk x — cos k y , and g(q) has to be negative, at least for q — > 0. The resulting Lan- 
dau function / kk / = /kk'(0) captures qualitatively the most pronounced features of the 
Landau function obtained from renormalization group calculations [3] and perturbation 
theory |TT] for the two-dimensional Hubbard model near van Hove filling. In particular, it 
is repulsive for momenta near two different van Hove points and attractive for momenta 
near a common one. 

We compute the two-particle vertex function T from / by summing the series of bubble 
chains sketched in Fig. 1. By virtue of the separable structure of the interaction, the series 
can be summed algebraically yielding 



where Hd is the particle-hole bubble with a form factor d^ at the vertices, that is 



Here / is the Fermi function. For small q and lo the d-wave bubble has the following 
asymptotic behavior: ILj(q, 0) = — N[ + a(q) q 2 + 0(q 4 ) and ImILj(q, u) — > —c(q)uj/q 
for q, to —> and cu/q — > 0, where iVj > is a weighted density of states, with each state 
at the Fermi level weighted by the squared form factor c? k , while a(q) and c(q) are real 
coefficients depending only on the direction q of q, but not on its length. The coefficient 
c(q) is always positive, but the sign of a(q) depends on q and on the choice of ek, d k and 
/i. 



r kk '(q,^) 



i-#(q) n d (q,u;) 



(4) 
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The PI sets in when the denominator of Eq. (4) vanishes at zero frequency and van- 
ishing momenta, that is lim q ^o <?(q) Il^q, 0) = —g(0)N[ = 1 at the critical point. We 
assume that the q-dependence of g(q) is such that g(q) ILi(q, 0) < 1 for q ^ at that 
point. Otherwise an instability with a finite q- vector would set in first. For a small 
q-vector this would lead to a phase where the symmetry-breaking Fermi surface deforma- 
tions are slowly modulated across the system. We leave this case for future studies. 

At the critical point the two-particle vertex has the following asymptotic form for 
small q, uj and small oo/q: 



with a(q) = a(q) — N[ g"(q)/g(0), where g"(q) is the second derivative of g(q) with 
respect to q for q — > 0. The denominator of the vertex has the same form as for other 
familiar critical points in metals, namely at the boundary to ferromagnetic ^21 EI] or phase 
separated ^3] states. Peculiar to the d-wave PI are the form factors in the numerator. 

We stress that the PI obeys the symmetry conditions for a continuous phase transition, 
since odd powers of the order parameter are excluded by symmetry from the correspond- 
ing Landau theory, in contrast to the case of phase separation, where cubic terms usually 
drive a first order transition (except for special cases where the cubic term is tuned to 
zero). Note also that the Pomeranchuk instability is not inhibited by long-range Coulomb 
forces, since volume conserving Fermi surface deformations do not generate charge inho- 
mogeneities. 

To estimate the decay rate for single-particle excitations in the presence of critical 
Fermi surface fluctuations near the quantum critical point, we have computed the imagi- 
nary part of the self-energy to first order in r k k' (q, a;) , evaluating the Feynman diagram 
in Fig. 2. For momenta on the Fermi surface, the result is 



r kk /(q,a;) 



<i k d k > 



(6) 



ic(q) uj jq — ct(q) q 2 



ImE(k F , e) oc dl F e 2/3 



(7) 
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for small energies e at T = and 

Im£(k F ,0) oc dl F T 2/3 (8) 

at low finite temperatures. Note that the integral leading to the above result is dominated 
by small momentum transfers q and uj <C Vk F q, thus justifying the asymptotic expansion 
of the vertex function. Furthermore, the dominant contributions to ImS at the Fermi 
vector kp come from (small) q-vectors that are tangential to the Fermi surface at kp. 

Not unexpectedly, the decay rate has the same energy dependence as for the quantum 
critical point near phase separation in two dimensions. Jl] Different is however the d-wave 
form factor making the decay rate strongly anisotropic. The decay rate is strongest near 
the van Hove points, while the leading terms vanish on the diagonal of the Brillouin zone. 
Subleading terms will produce at least conventional Fermi liquid decay rates (T 2 log T) on 
the diagonal, but faster decay (intermediate between T 2 / 3 and T 2 log T) may be obtained 
due to higher order processes and interactions with large momentum transfers, which 
couple different parts of the Fermi surface. 

The singular self-energy in the critical region will drastically modify the one-particle 
propagator at low energies. However, this does not invalidate the calculation of the two- 
particle vertex from bubbles with bare propagators, since singular self-energy corrections 
generated by strong forward scattering are cancelled by corresponding vertex corrections 
in the polarization bubble. P5J A more subtle point is whether higher order corrections 
to the self-energy contribution in Fig. 2 modify the power-laws (7) and (8) for the decay 
rate. The same question has been discussed at length for fermions coupled to a gauge 
field. [15] where the lowest order calculation also yields a power-law with exponent 2/3.[T6] 
A detailed analysis by Altshuler et al. [T2j indicated that this result is not changed by 
other terms. A recent renormalization group analysis further supported the validity of 
the leading power-law for the gauge theory, and also for the critical point near phase 
separation in two dimensions. [IHj Even if the above power-law would be modified by 
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higher order corrections, it is clear that Fermi liquid behavior cannot be restored, and the 
decay rate will remain large and anisotropic in any case. 

A strongly anisotropic decay rate following a power-law with exponent 2/3 has recently 
been derived for an isotropic continuum (not lattice) version of model (2). [7] That result 
was obtained for the symmetry-broken " nematic" phase, and the anisotropy of the decay 
rate arises from the anisotropy of the symmetry-broken state and its collective modes. In 
the symmetric phase of an isotropic liquid the decay rates are of course isotropic. At the 
quantum critical point the decay rate of the isotropic liquid also obeys a power law with 
exponent 2/3, but now isotropically over the whole Fermi surface. 

Could soft Fermi surfaces and critical Fermi surface fluctuations play a role in cuprate 
superconductors? Due to the coupling of electron and lattice degrees of freedom a 
symmetry-breaking Fermi surface deformation is generally accompanied by a lattice dis- 
tortion, and vice versa. Structural transitions which reduce the lattice symmetry of the 
cuprate-planes are quite frequent in cuprates. Close to a PI of the electronic system, elec- 
tronic properties can be expected to react unusually strongly to slight lattice distortions. 
Such "overreactions" of electronic properties have indeed been observed in several cuprate 
compounds. JH] I n particular, a slight orthorhombicity of the lattice structure would lead 
to a relatively strong orthorhombic distortion of the Fermi surface. Yamase and Kohno 
[2Hj invoked this idea to explain peculiarities of magnetic excitations in cuprates. 

Large Fermi surface fluctuations could be at least partially responsible for the non- 
Fermi liquid behavior observed in cuprates at optimal doping. In our model calculation 
we have obtained a strongly anisotropic anomalously large decay rate for single-particle 
excitations. Large anisotropic decay rates have been frequently inferred from the linewidth 
observed in photoemission experiments on optimally doped cuprates. [21] However, recent 
experiments with higher resolution revealed that bilayer splitting of the bands has to be 
taken into account in the data analysis, such that the intrinsic linewidth near the van Hove 
points may be considerably smaller and less anisotropic than previously expected. (22] 
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As to the temperature dependence of the decay rate, it is presently hard to reliably 
discriminate a T-linear from a T 2 / 3 behavior of the intrinsic linewidth extracted from the 
experimental data. Concerning transport, an anisotropic scattering rate with nodes on 
the diagonal can very naturally account for the pronounced anisotropy between intra- 
and inter-plane mobility of charge carriers, as pointed out by Ioffe and Millis (23] in their 
phenomenological "cold spot" scenario. According to their idea, the intra-plane transport 
is dominated by quasi-particles with a long life-time near the diagonal of the Brillouin 
zone, while these carriers are not available for inter-plane transport, since transverse 
hopping amplitudes vanish on the diagonal. To analyze transport properties near a PI 
one has to face the notorious difficulties of transport theory in a quantum critical regime. 
Even for rough estimates of temperature dependences within a semiclassical Boltzmann 
approach one would have to compute subleading corrections to the decay rates, since 
those limit the life-time of excitations on the Brillouin zone diagonal, and thus the in- 
plane conductivities. Furthermore, regular interactions or interactions peaked at finite 
q- vectors (leading to "hot spots" on the Fermi surface) should also play a role, and may 
lead to an interesting but complicated interplay with Fermi surface fluctuations. 

Since the PI breaks the orientational symmetry of the lattice, it is natural to consider 
a possible connection with the tendency towards stripe formation, which has been ex- 
tensively discussed in the context of cuprate superconductors. [25 Stripes also break the 
translation invariance in addition to orientational symmetry, and their formation requires 
interactions with large momentum transfers, such as antiferromagnetic interactions. The 
possibility of direct transitions between states with broken orientational order and stripe 
states, which has been envisaged already by Kivelson et al.jH], is an interesting subject 
for future studies. 

In summary, an electron system close to a PI is characterized by a soft Fermi sur- 
face, which can react strongly to a slight change of the lattice structure, and exhibits 
strong collective Fermi surface fluctuations. In two-dimensional systems these fluctu- 
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ations lead to large anisotropic decay rates for single-particle excitations and thus to a 
breakdown of Fermi liquid theory. It will be interesting to further explore the consequences 
of Pomeranchuk criticality, especially for charge and heat transport, for magnetism, and 
for superconductivity. 

Acknowledgements: 

We are very grateful to Claudio Castellani for contributing several valuable ideas and in- 
sights, especially at the initial stage of this work. We also thank B. Buchner, B. Keimer, C. 
Di Castro, G. Khaliullin, G. Kotliar, V. Oganesyan, M. Salmhofer, F. Wegner, H. Yamase, 
and R. Zeyher for useful discussions. This work was initiated at the workshop "Realistic 
theories of correlated electron materials" in fall 2002 at ITP-UCSB, and supported by the 
DFG-grant Me 1255/6-1. 

References 

[1] I.J. Pomeranchuk, Sov. Phys. JETP 8, 361 (1958). 

[2] H. Yamase and H. Kohno, J. Phys. Soc. Jpn. 69, 332 (2000); 69, 2151 (2000). 

[3] C.J. Halboth and W. Metzner, Phys. Rev. Lett. 85, 5162 (2000). 

[4] I. Grote, E. Kdrding, and F. Wegner, J. Low Temp. Phys. 126, 1385 (2002); V. 
Hankevych, I. Grote, and F. Wegner, Phys. Rev. B 66, 094516 (2002). 

[5] A. Neumayr and W. Metzner, Phys. Rev. B 67, 035112 (2003). 

[6] B. Valenzuela and M.A.H. Vozmediano, Phys. Rev. B 63 153103 (2001). 

[7] V. Oganesyan, S.A. Kivelson, and E. Fradkin, Phys. Rev. B 64, 195109 (2001). 

[8] S.A. Kivelson, E. Fradkin, and V.J. Emery, Nature 393, 550 (1998). 



9 



[9] Although the discontinuity of the momentum distribution function is smeared by the 
gap in a superconductor, a Fermi surface can still be defined as the locus of zeros of 
the diagonal components of the inverse Nambu propagator in momentum space. 

[10] M. Vojta, Y. Zhang, and S. Sachdev, Phys. Rev. Lett. 85, 4940 (2000); Int. J. Mod. 
Phys. B 14, 3719 (2000). 

[11] PA. Frigeri, C. Honerkamp, and T.M. Rice, Eur. Phys. J. B 28, 61 (2002). 

[12] J.A. Hertz, Phys. Rev. B 14, 1165 (1976). 

[13] A.J. Millis, Phys. Rev. B 48, 7183 (1993). 

[14] C. Castellani, C. Di Castro, and M. Grilli, Phys. Rev. Lett. 75, 4650 (1995). 

[15] For a review on Fermi systems with strong forward scattering, see W. Metzner, C. 
Castellani, and C. Di Castro, Adv. Phys. 47, 317 (1998). 

[16] PA. Lee, Phys. Rev. Lett. 63, 680 (1989); B. Blok and H. Monien, Phys. Rev. B 47, 
3454 (1993). 

[17] B.L. Altshuler, L.B. Ioffe, and A.J. Millis, Phys. Rev. B 50, 14048 (1994). 

[18] C. Castellani, S. Caprara, C. Di Castro, and A. Maccarone, Nucl. Phys. B 594, 747 
(2001). 

[19] J.D. Axe et al, Phys. Rev. Lett. 62, 2751 (1989); B. Biichner et al., Phys. Rev. Lett. 
73, 1841 (1994). 

[20] H. Yamase and H. Kohno, J. Phys. Soc. Jpn. 70, 2733 (2001); ibid 71, 1154 (2002). 

[21] See, for example, T. Valla et al., Phys. Rev. Lett. 85, 828 (2000). 

[22] For a recent review on photoemission in cuprates, see A. Damascelli, Z. Hussain, and 
Z.X. Shen, Rev. Mod. Phys. 75, 473 (2003). 

10 



[23] L.B. Ioffe, and A.J. Millis, Phys. Rev. B 58, 11631 (1998). 

[24] See, for example, Proc. of 3rd Int. Conf. on Stripes and High T c superconductivity, 
Int. J. Mod. Phys. B 14, Nos. 29-31 (2000). 



11 



f f 




+ 



Figure 1: Series of bubble chains contributing to the two-particle vertex T. 
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Figure 2: Feynman diagram relating the self-energy E to the two-particle vertex T. 
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